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This paper describes high-energy neutron inelastic scattering measurements of propagating mag- 
netic excitations in PrBa2Cu306+:c (x m 0.2 and 0.93). The measurements probe the acoustic and 
optic modes of the antiferromagnetically-ordered copper-oxygen bilayers in the energy range 50- 
150 meV. The observed magnon dispersion can be described satisfactorily in this energy range by a 
spin wave model including intra- and inter-layer nearest-neighbour exchange constants Jii and J±. 
We find Jy = 127 ± lOmeV and J± = 5.5 ± 0.9 meV. The value of Jy is virtually the same as that 
found in YBa2Cu3C>6.2, but J± is a factor of two smaller. To within experimental error the values 
of Jy and J± for Prl^CusOe+i do not vary with oxygen doping. 

PACS numbers: PACS numbers: 74.72.Bk, 75.30.Ds, 75.50.Ee, 78.70.Nx 



I. INTRODUCTION 

One of the unusual characteristics of the cuprate su- 
perconductors is the relative insensitivity of the super- 
conducting properties to the presence of magnetic rare- 
earth ions. In the -ffi^CusOe+z family (R = Y or 
rare earth) superconductivity occurs at temperatures as 
high as 95 K, and magnetic ordering of the R sublat- 
tice generally coexists with superconductivity at temper- 
atures around 2 K or below i Antiferromagnetic (AFM) 
ordering of the bilayer Cu spins is observed in non- 
superconducting samples with low oxygen doping levels 
(typically x < 0.4), but superconducting samples do not 
exhibit any conventional form of Cu magnetic order p* 

A striking exception to this norm is the case of 
PrBa2Cu306+a; (hereafter PrBC06+x), which exhibits 
anomalous electrical and magnetic properties in compar- 
ison to other i?BCO compounds. 3,4 Superconductivity is 
not found for any x in samples of PrBC06+x prepared 
by standard methods^, and the transition temperature 
Tn ~ 300 K for AF order of the Cu spins depends only 
weakly on xfi Magnetic ordering of the Pr sublattice 
takes place below a temperature Tp r varying from 11 K 
(x = 0) to 17K (x = much higher than the rare 

earth magnetic ordering temperatures in superconduct- 
ing i?BCO compounds. Unusual magnetic structures are 
observed below Tp r due to magnetic coupling between 
the Cu and Pr sublatticesiiS 

Many of these anomalous features of PrBC06+a; re- 
main a puzzle. Models for the electronic structure have 
indicated a tendency for Pr 4/ - O 2p hybridization to 
cause a localization of doped holes in O 2p v orbitals, thus 
inhibiting superconductivity. 11 Such hybridization would 



also be expected to influence the magnetic couplings in 
PrBC06+i£, and so measurements of the exchange inter- 
actions could yield information on the underlying elec- 
tronic structure. 

Over a number of years we have undertaken a sys- 
tematic investigation of the magnetic excitations of 
PrBC06+a; by neutron inelastic scatteringiiSii 3 - 1 ^ The 
aim has been to identify the important magnetic cou- 
plings and to see how they vary with doping. Broadly 
speaking, the magnetic excitations fall into two cat- 
egories: (i) transitions between levels of the Pr 4/ 
electrons split by the local crystalline electric field 
(CEF)?i2*iii and (ii) Cu spin wave excitations. There are 
also non-trivial effects on the excitations due to the mag- 
netic coupling between the Pr and Cu sub-systems^ 1 ^ 

In our most recent work we have investigated the 
wavevector dependence of the magnetic excitations us- 
ing a single crystal sample of PrBC06+a; prepared first 
in an oxygen-deficient state (x « 0.20) and later treated 
to produce an 'opti mally -doped' oxygen content {x « 
0.93). Refs. E3 and llj describe measurements at low 
(< lOmeV) and intermediate (45-65 meV) energies, and 
present an analysis of the spectra in terms of a spin 
model for the coupled Cu-Pr system. Here we report 
measurements of the Cu spin excitations at higher en- 
ergies (50-150 meV) by time-of-flight neutron scattering. 
We determine the principal intra- and inter-layer Cu- 
Cu exchange interactions Ju and Jj_ in PrBC06+:r and 
make a comparison with the corresponding parameters in 
YBCOG+rr. Perhaps the most interesting new finding is 
that Jii and Jj_ in PrBC06+x are independent of oxygen 
doping. 
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II. EXPERIMENTAL DETAILS 

The experiments were performed on the same single 
crystal of PrBC06+a: as used for the measurements de- 
scribed in Refs. lip an d 113. and also for the phonon study 
reported in Ref. [jja- The as-grown crystal had a mass 
of 2 g, and a mosaic spread of approximately 1 deg (full 
width at half maximum) measured by neutron diffrac- 
tion. The crystal was subjected to two annealing treat- 
ments to control the oxygen content. First, the as-grown 
crystal was reduced at 700°C in a flow of 99.998% ar- 
gon for 100 h and quenched to room temperature. After 
the measurements on this 'under-doped' crystal had been 
completed the crystal was then annealed in pure oxygen 
over a period of 60 days while cooling in steps from 600°C 
to 450° C, with progressively longer dwell times at each 
step as the temperature decreased. After the second an- 
neal the mosaic of the crystal was found to have increased 
to approximately 3 deg. 

According to studies of oxygen content in YBC06+X 
as a function of annealing conditions (oxygen partial 
pressure and temperature) the two treatments applied 
here yield oxygen contents of x < 0.2 and x « 0.93. 17 
Although the oxygen ordering properties of YBC06+X 
and PrBC06+x are qualitatively different their oxygen 
equilibrium isobar coincides at 0.21 atm up to 650°C. 18 
Therefore it is realistic to assume that this is also the 
case for other oxygen partial pressures. However, as 
a further check of the oxygen content we measured 
the magnetic ordering behaviour of the crystal after 
each anneal by neutron diffraction. Studies carried out 
on polycrystalline PrBC06+a: have shown by various 
tcchniquesSiLSiS that the Pr magnetic ordering temper- 
ature increases from Tp T — 1 1 K for vacuum-annealed 
samples (x « 0), to Tp r = 17 K for 'optimally-doped' 
material (x « 0.93). The results for our crystal are 
illustrated in Fig. ^ which shows the temperature de- 
pendence of the (0.5,0.5,0) Bragg reflection whose ap- 
pearance is known to signal the onset of Pr magnetic 
order piU-!^-^ The Pr ordering temperatures are deter- 
mined to be Tp r = 13 ± 0.5 K after reduction, and 
Tp r = 17.5 ± 0.5 K after oxygen annealing. These values 
are consistent with the literature results for the expected 
oxygen contents x w 0.2 and x ~ 0.93, respectively, and 
we will henceforth refer to the crystals as PrBC06.2 and 
PrBC06.93. The sharpness of the transitions indicates 
a homogeneous oxygen distribution. The (0.5,0.5,0) in- 
tensity curve for PrBC06.93 also shows an anomaly at 
T% = 12.5 ± 0.5. This is consistent with a spin transi- 
tion first observed by Uma et a/.fSi which corresponds to 
a change in the c axis stacking sequence of the coupled 
Pr-Cu magnetic structured 

The neutron inelastic scattering measurements re- 
ported here were performed on the MAPS (PrBC06.2) 
and HET (PrBC06.93) spectrometers at the ISIS spalla- 
tion neutron source. MAPS and HET are time-of-flight 
spectrometers equipped with pixellated area detectors 
situated 6 m (MAPS) and 4 m (HET) from the sample po- 
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FIG. 1: (Color online). Temperature dependence of the 
(0.5, 0.5, 0) magnetic Bragg peak of PrBC06+a;. This peak 
signals magnetic order on the Pr sublattice. The curves were 
recorded after the annealing treatments described in the text, 
which resulted in oxygen contents of x ~ 0.2 and x ps 0.93. 
The measurements were made by neutron diffraction on the 
triple-axis spectrometers TAS6 at Ris0 National Laboratory 
and IN14 at the Institut Laue-Langevin. The maximum in- 
tensity has been scaled to unity for the sake of comparison. 



sition, and an incident beam divergence of approximately 
0.5°. These design features ensure very good wavevector 
(Q) resolution, an essential experimental requirement in 
this work because of the very steep Cu spin wave disper- 
sion. 

The crystal was aligned with the c axis approximately 
parallel to the incident neutron beam direction. This 
arrangement means (i) that the surfaces in Q space cor- 
responding to constant neutron energy transfer are ap- 
proximately parallel to the a*b* plane in the reciprocal 
lattice of the crystal, and (ii) that the component of Q 
parallel to c* varies with energy transfer. The natural 
way to visualize the data is then to integrate the signal 
recorded in each pixel of the detector bank over a band 
of energies and project the results on to the a*b* plane. 

Measurements were performed with fixed incident neu- 
tron energies of 80 meV, 200 meV and 300 meV on MAPS, 
and 200 me V on HET. Most of the data were collected 
with the crystal at a temperature of 20 K, but some runs 
were repeated at 300 K. The runs on MAPS were of 2-3 
days duration, and the runs on HET were of 5 days dura- 
tion for each of 2 settings of the crystal. These run times 
are given for an average proton current of 160 fiA. In prin- 
ciple, the range of incident energies employed on MAPS 
allowed us to study the spin excitations up to ~ 260 meV, 
which is the anticipated maximum in the one-magnon 
spectrum (see later). However, owing to the relatively 
small size of the crystal the count rate was too low to 
obtain statistically meaningful data above 150 meV. The 
intensity was converted into an absolute cross section by 
comparison with the scattering from a standard sample 
of vanadium^i The presented spectra are the partial dif- 
ferential cross section d 2 a/dttdE{ per formula unit (f.u.) 
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FIG. 2: (Color online). The AF1 magnetic structure. The 
diagram depicts one unit cell of PrBC06+:r. Only the Cu 
atoms are shown, and ordered magnetic moments are carried 
only by the Cu atoms in the bilayer. The constants Jy and 
J± are the nearest-neighbour intra- and inter-layer exchange 
parameters, respectively. 



multiplied by the factor ki/kf, where fc; and k{ are the 
initial and final neutron wavevectors, and Ef is the final 
energyiaa 

III. CU SPIN EXCITATION SPECTRUM 

As mentioned above, the bilayer Cu spins in 
PrBC06+a; order antiferromagnetically at a temperature 
Tn near room temperature and slightly dependent on x. 
At temperatures below Tn but above Tp r the bilayer Cu 
spins are aligned antiparallel to their nearest neighbours 
along all three crystallographic directions, as shown in 
Fig. El There is no ordered moment on the Cu site in 
the Cu-0 chains. This is the well known AFI structure, 
the same as found in undcrdoped YBCO6+2;. Our in- 
terpretation of the data can therefore be guided by the 
spin wave model used to describe the spin excitations in 
AFI-ordered YBCQ6+a;i 23 i 24 i 25 

For the high-energy spin excitations it is a good ap- 
proximation to regard the Cu spin arrangement as a 
square-lattice bilayer antiferromagnet. This neglects the 
slight orthorhombic distortion (if present) and the weak 
coupling between bilayers in adjacent unit cells. The sim- 
plest spin Hamilatonian that describes this system is then 

h = Jwj^i ■ s 3 - + j±j2 s i ■ s 3 r , (!) 

(ij) W) 

The summations in Eq. Q are over pairs of nearest- 



neighbour Cu spins in the same layer (first term) and 
on adjacent layers (second term). Each pair of spins is 
counted only once. The constants Jm and Jj_ are the 
intra- and inter-layer exchange parameters, respectively. 

The magnon spectrum derived from Eq. JJl has two 
branches, differing according to whether the spins on 
adjacent layers in the bilayer rotate in the same sense 
(acoustic modes) about their average direction, or in the 
opposite sense (optic modes). The dispersion of these 
branches is given by2i2^ 

MQ) = 2J||{1 - 7 2 (Q) + (Jx/2J||)[1 ± 7(Q)]} 1/2 , (2) 

where + and — correspond to the acoustic and optic 
modes, respectively, and 

7(Q) = ^{cos(Q x a) + cos(Q y a)}. (3) 

Close to the 2D AFM zone centres, e.g. Q = (|, |) in 
units of 2n/a, the acoustic branch is linear with an initial 
slope proportional to J|| . The optic branch has an energy 
gap of 2y / (J|| Jj_) at the AFM zone centres. Because of 
the bilayer structure in the unit cell (Fig. |2J) the cross 
section for scattering from acoustic and optic magnons 
contains factors of sin 2 (;zQzc/2) and cos 2 (zQ z c/2), re- 
spectively, where z is the inter-layer spacing as a fraction 
of the c lattice parameter (z — 0.295 for PrBC06+x). 



IV. RESULTS 

Figures|3]andE]present some examples of neutron scat- 
tering data collected on the MAPS and HET spectrom- 
eters, respectively. Figs. Ola) and (b) are maps of the 
intensity from PrBC06.2 averaged over the energy range 
85-105 meV and 130-160 meV, respectively, projected on 
the (h, k) plane. For these energies we could not detect 
any difference between the signal at 20 K and at 300 K, 
so data collected at these two temperatures were aver- 
aged to improve statistics. The images focus on the re- 
gion of 2D reciprocal space in the vicinity of the AFM 
zone centre (0.5,0.5), and contain intensity in excess of 
background distributed around (0.5,0.5). This signal, 
which expands away from (0.5,0.5) with increasing en- 
ergy, is consistent with the scattering from AFM spin 
waves. According to Eq.[2 the scattering from spin wave 
excitations of the bilayer Cu spins is expected to be vir- 
tually isotropic in the plane (i.e. a ring of scattering) be- 
cause at these energies the spin wave wavevector extends 
only about halfway to the Brillouin zone boundary. The 
experimental signal is consistent with a ring centred on 
(0.5,0.5) once the statistical scatter of the data is taken 
into consideration. 

Under the conditions of the measurement the out-of- 
plane wavevector component at E — 95meV is Q z = 5.3 
(in units of 2ir/c), and at E = 145 meV it is Q z = 6.6. 
Maxima in the bilayer structure factors for acoustic and 
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FIG. 3: (Color online). Neutron scattering from 
PrBa2Cu306+:r (x « 0.2) measured on the MAPS spectrome- 
ter, (a) and (b) are constant-energy slices in which the inten- 
sity has been averaged over a band of energies from (a) 85 meV 
to 105 meV, and (b) 130 meV to 160 meV. (c) and (d) are cuts 
taken from the slices in (a) and (b) parallel to (h, 0) passing 
through (0.5,0.5). The solid lines are the results of fitting 
the data to two Gaussian functions on a sloping background. 
The incident neutron energies used in the measurements were 
200 meV [(a) and (c)] and 300 meV [(b) and (d)]. Data from 
the four equivalent Brillouin zones (±0.5, ±0.5) have been av- 
eraged. 
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FIG. 4: (Color online). Neutron scattering from 
PrBa2Cu306+i (x w 0.93) measured on the HET spectrom- 
eter with an incident energy of 200 meV. (a) is a constant- 
energy slice in which the intensity has been averaged over 
the range 90-110 meV, and (b) is a cut taken parallel to the 
(—h,h) direction passing through (0.5,0.5). The solid lines 
are the results of fitting the data to two Gaussian functions 
on a sloping background. 



optic magnons are found at Q z — 5.1 and 6.8, respec- 
tively, so the scattering is mainly from acoustic modes in 
Fig. 01a) and from optic modes in Fig. E^b). 

Figs. |3Jc) and (d) show linear scans extracted from the 
intensity maps in Figs. Ofa) and (b). The scans are made 
parallel to (h, 0) passing through the AFM zone centre 
at (0.5,0.5), and are averaged over a range Afc = 0.08 
(95meV) and Ak = 0.1 (145 meV) in the (0, k) direc- 
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FIG. 5: (Color online). In-plane magnon dispersion of 
PrBa2Cu306+a;. Circles and squares denote measurements 
of the acoustic and optic branches, respectively. Open and 
closed symbols are for x — 0.2 and x = 0.93, respectively. 
The data points for the optic mode gap at h = 0.5 for x — 0.2 
and x — 0.93 are virtually coincident, as indicated. Vertical 
'error bars' indicate the energy range over which the data were 
averaged, whereas horizontal error bars show the experimen- 
tal uncertainty in the radius of the dispersion surface. The 
lines are the spin wave dispersion relations for a bilayer anti- 
ferromagnet calculated from Eqs.|5] and |3] with Jm = 127 meV 
and J± — 5.5 meV. 



tion. Each scan contains two peaks, one for each inter- 
section of the scan with the spin wave dispersion surface 
which we assume to be a ring. By fitting a pair of Gaus- 
sian functions on a sloping background to the peaks in 
these and similar scans in other directions and from other 
runs, allowing for the curvature of the dispersion surface 
over the averaging width Ak, we found the radius of the 
ring as a function of energy, and hence arrived at the 
in-plane spin wave dispersion shown in Fig. [3] We have 
chosen to plot the dispersion along the (h, h, 0) direction 
for consistency with published data on YBCO^Si^ but 
we reiterate that the spin wave dispersion is expected 
to be virtually isotropic within the range of wavevectors 
probed in our measurement. The points corresponding 
to mainly acoustic modes and mainly optic modes are 
shown with different symbols on Fig. [21 and we have also 
included the optic mode gap at the AFM zone centre 
(53 ± 2 meV) determined in a previous measurement on 
the same PrBC06.2 crystal^ 

The magnon dispersion data were compared with the 
curves for the acoustic and optic spin wave branches cal- 
culated from Eqs.[21and[31 The two exchange parameters 
were adjusted to give the best overall agreement with the 
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experimental data, and after consideration of the experi- 
mental uncertainties we obtained Ju = 127 ± lOmeV and 
J± = 5.5 ± 0.9 meV. The best-fit dispersion curves are 
plotted on Fig. [SJ 

The data collected on oxygenated PrBC06.93 were 
more limited than those just described for PrBC06.2. 
This is because the detector area on HET is consider- 
ably smaller than on MAPS, which restricted the mea- 
surement to one Brillouin zone and required two slightly 
different settings of the crystal to cover enough of the 
zone to include all the signal. Nevertheless, the constant- 
energy slice shown in Fig. 0Ja) has Q z — 5.3, and so 
corresponds to an almost pure acoustic mode. Fig. 0Jb) 
shows a cut through the antiferromagnetic point (0.5, 0.5) 
parallel to (— h, h). This direction was chosen for the cut 
because the data recorded in the detector extend fur- 
thest along this diagonal, as can be seen in Fig. 0Jb). 
From the two-Gaussian fit we obtained the point on the 
acoustic spin wave dispersion curve shown on Fig. |SJ To 
within experimental error (~ 10%) the PrBC06.93 da- 
tum lies on the PrBC06.2 dispersion curve. In an ear- 
lier experiment^* we determined the optic mode gap for 
PrBC06.93 and found it to be 54 ± lmeV, again consis- 
tent with that for PrBC06.2. 



V. DISCUSSION 

At the start of this work we set out to determine the 
Cu-Cu exchange parameters Jii and J± for PrBC06+ir 
and to compare them with the corresponding parameters 
for YBCO6+2;. The result Jy = 127 ± lOmeV for the 
intra-layer exchange parameter of PrBC06.2 found here 
is the same to within experimental error as the values 
J|| = 125 ±5 me V 25 and Jy = 1 20 ±20 meV^ determined 
by neutron inelastic scattering for YBC06.15. There is 
also good agreement with the values J11 « 115meV for 
PrBCOG+a; {x w 0), Refjp and Jy m 120 meV for 
YBC06+a; (x ~ 0), Ref. |23, derived from two-magnon 
Raman scattering, which provides some support for the 
models used to describe the rather broad two-magnon 
Raman spectra. 

On the other hand, the values of the inter-layer ex- 
change parameter for PrBC06+ir and YBC06+X are 
significantly different. We obtain J± = 5.5 ± 0.9 meV 
for PrBC06.2, compared with J± = 9 — lOmeV for 
YBC06.2 (Ref.lH and J ± = 11 ± 2meV for YBC06.15 
(Ref. |25T> . Our own measurements on a crystal of 
YBC06.2 gave J± = 13 ± 2meV. One factor that might 
influence J± is the inter-layer separation, which is about 
4% larger in PrBC06+a; than in YBC06+X reflecting 
the larger size of the Pr 3+ ion relative to Y 3+ . It would 
be surprising, however, if this accounts for a factor 2 dif- 
ference in J± between PrBCO and YBCO. 

As mentioned in the introduction, it has been pro- 
posed that superconductivity is suppressed in PrBCO be- 
cause holes become localized in hybridized Pr-0 bondsii. 
Naively, one might expect such a change in the electronic 



structure to have an influence on the Cu-Cu exchange in- 
teractions. The invariance of Jy to the replacement of Y 
by Pr implies that the in-plane AFM superexchange in- 
teraction is unaffected by any changes in the electronic 
structure associated with Pr 4/ - O 2p hydridization. 
This does not present any obvious difficulties since the 
superexchange is mediated by the O 2p a orbitals which 
lie in the layers, whereas the proposed hybridization 
schemeAi involves O 2p^ orbitals which are oriented per- 
pendicular to the layers. Following the same reasoning, 
it is tempting to connect the large difference in J± be- 
tween PrBCO and YBCO as evidence for hybridization- 
induced changes in the electronic structure within the bi- 
layer. However, if the J± exchange were mediated by hy- 
bridized Pr-0 bonds that also accommodate doped holes 
then it is difficult to explain the apparent insensitivity of 
J± to doping. A fuller analysis of these observations will 
require a proper understanding of the mechanism for the 
inter-layer exchange in the bilayer cuprates. What we 
can say, though, is that the substitution of Y by Pr does 
have a very significant influence on the inter-layer cou- 
pling. 

Before concluding it is worth mentioning that although 
we have found no difference between the spin wave 
spectrum of PrBC06.2 and PrBC06.93 below 100 meV, 
there may still be differences at higher energies. In- 
deed, the two-magnon Raman peak, which is sensitive to 
short wavevelength spin fluctuations, was found shifted 
to lower energies in spectra from PrBC07 compared 
with PrBC06f2LSi leading to a significantly smaller value 
Jii i=a 95 meV than obtained from our neutron scattering 
measurements. One explanation for this apparent dis- 
crepancy could be that the magnetic excitations in hole- 
doped PrBCO may be influenced by magnon-hole inter- 
actions at higher energies. Neutron scattering measure- 
ments of the magnetic excitations up to the one-magnon 
zone boundary could therefore be very interesting. 



VI. CONCLUSIONS 



In this work we have shown that the Cu spin excitation 
spectrum in PrBCO is well described up to ^150 me V by 
the spin wave model for a bilayer antiferromagnet. In this 
energy range the spin excitations do not change with dop- 
ing. The in-plane superexchange parameter Jy has the 
same value in undoped PrBCO as in undoped YBCO, 
whereas the inter-layer exchange parameter J± is a factor 
of 2 smaller in PrBCO than in YBCO. An understanding 
of these results could provide useful insight into the dif- 
ferences between the electronic structure of PrBCO and 
YBCO, and hence into the question of why YBCO has a 
superconducting ground state whereas PrBCO does not. 
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